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Abstract 

We establish rigorous upper and lower bounds for the speed of pulled fronts with a cutoff. We 
show that the Brunet-Derrida formula corresponds to the leading order expansion in the cut-off 
parameter of both the upper and lower bounds. For sufficiently large cut-off parameter the Brunet- 
Derrida formula lies outside the allowed band determined from the bounds. If nonlinearities are 
neglected the upper and lower bounds coincide and are the exact linear speed for all values of the 
cut-off parameter. 
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The reaction diffusion equation 



iH = u xx + f(u) (f ) 

provides a simple description of phenomena in fields such as population dynamics, chemical 



> population a 1 



It is one of the 



reactions, flame propagation, fluids, QCD, among others 
simplest models which shows how a small perturbation to an unstable state develops into a 
moving front joining a stable to an unstable state. The reaction term f(u) satisfies different 
conditions depending on the physical problem of interest. One of the first, and most studied 
cases, is the Fisher reaction term f(u) = u(l — u) for which the asymptotic speed of the 
propagating front is c = 2, a value determined from linear considerations. A more general 
case was studied by Kolmogorov, Petrovskii and Piscounov (KPP)^ who showed that for 
all reaction terms which satisfy the KPP condition 

f(u) > 0, /(0) = /(l) = 0, /(«) < f(0)u (2) 

the asymptotic speed of the front joining the stable u — 1 point to the unstable u = point 
is given by 



CKPP = 2 y/f (0). 

These fronts are called pulled since it is the leading edge of the front which determines the 
velocity of propagation. In the rest of this work we assume that /'(0) = 1. The evolution of 
localized initial conditions for general reaction terms, and rigorous properties of the fronts 



were studied by Aronson and Weinberger j^]. The asymptotic speed of the front for all 
reaction terms can be found from the integral variational principle js] 

c 2 = sup 2 f f^fff" (3) 



where the supremum is taken over all positive monotonic decreasing functions g(u) for which 
the integrals exist and where h(u) = —g'(u). The supremum is always attained for reaction 
terms which are not pulled. 

Two effects not included in the classical reaction diffusion equation (pQ), are the effect of 
noise and the effect of a finite number N of diffusive particles. It was shown by Brunet and 
Derrida that such effects can be simulated by introducing a cut-off in the reaction term. In 
the case of noise the cut-off parameter measures the amplitude of the noise while in the case 
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of finite number of N diffusing particles the cut-off parameter e = 1/N. There is substantial 
numerical evidence that introducing a cut-off in the reaction terms reproduces accurately 
the effect of noise and finiteness in the number of diffusing particles js]. 

By means of an asymptotic matching Brunet and Derrida showed that for a reaction term 
f(u) = it(f — it 2 ) a small cut-off changes the speed of the front to 



(hie) 2 ' 

In recent work it has been show that the Brunet-Derrida formula for the speed is correct to 

r~i 

0((\ne)~ 3 ) for a wider class of reaction terms [101 ]. 

The purpose of this work is to show that for reaction terms of the form f(u)Q(u — e) 
where / satisfies the KPP condition Eq. (J2J) and G is the step function, the speed c of the 
front with the cutoff satisfies 

2 sin(0*) - AO*) < c < 2 sin(0*), (5) 

with 

0* tan(0*) = -| ln(e)| (6) 

We see that for < e < 1, < 0* < 7r/2. The function A (</>*) depends on the nonlinear 
terms of the reaction function. For small e the series expansion of the upper bound cjjp is 

c UP = 2 sin(^) = 2 - + 0((ln e)~ 3 ) (7) 

The contribution of the nonlinearities, contained in the term A(0*), appears at 0((lne)~ 3 ), 
so that the leading order terms in the expansion of the upper and lower bounds give the 
Brunet-Derrida formula. If nonlinearities are neglected the value 2sin(0*) is the analog of 
the KPP value c = 2 for reaction terms which satisfy the KPP condition, but with a cutoff. 



In what follows we derive the bounds and apply them to the Fisher reaction term 11] 
f(u) = u — u 2 and to the reaction term studied by Brunet and Derrida f(u) = u — u 3 . The 
main tool to obtain the bounds is the variational principle for the speed. 

As shown in previous work [121 ]. we may perform the change variables u = u(s) where 
s = 1/g in Eq.(j3]) and write the variational expression for the speed as 



c 2 = su P 2 v 77 u ' JU v-v""/— ; ( 8) 
u(s) Jq° (du/ds) ds 
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where So — l/5 f (^ — 1) is an arbitrary parameter, 

F(u) = / f(q)dq. 
Jo 

and the supremum is taken over positive increasing functions u(s) such that u(0) = 0, 
u(sq) = 1 and for which all the integrals in (jSj) are finite. Therefore, for any suitable trial 
function u(s) we know that 

c2>2 F(l)/s + f°°F(u(s))/s*ds 

J°°(du/ds) 2 ds ' 1 ' 

Consider now reaction terms f(u) with a cut-off e of the form 

[ if < u < e 

m = - - 

I u - N[u) if e < m < 1, 
where N(u), the nonlinearity, is such that N(0) = N'(0) = 0. We find 



F(u) 



if < u < e 

u 2 /2 - e 2 /2 + F n (u) if e < u < 1, 



where F n (u) = — /" N{u)du. 

Assume now that /(m) satisfies the KPP criterion Eq.pj). Since f(u) < u, it follows that 
F(u) < where 

% if < u < e 

w 2 /2-e 2 /2 if e< u< 1, 



and therefore 



c<g[nj = sup2 — 2 . (10) 

«(s) J {du/ds) ds 

One can prove (rigorous details will be given elsewhere) that Q is bounded above and that 
there exists a function u(s) for which the supremum is attained. This function is the mono- 
tonic increasing solution to the Euler-Lagrange equation for Q satisfying the boundary con- 
ditions u(0) = 0,u(sq) = 1. One can also prove that the variational parameter sq is finite 
and u'(sq) = 0. In summary, the maximizing function for Q is the solution of 



subject to the boundary conditions 

u(0) = 0, u(s Q ) = 1 u'(s ) = u'(s) > 0, 

with the function and its derivative continuous at u = e. 
The solution to this problem is given by 



uis) 



s if < s < e 

cos(0(s)) if e < s < Sq, 



with 



A = y/e sec (0*), s = 1/e, 0(s) = - cot(0*) ln(s/e) — 0* (11) 



where 0* is the first positive solution of 

0* tan0* = ^| lne|. (12) 

The maximum of Q = Q[u] can be calculated easily. We obtain after performing the 
integrals, 

c 2 < Q[u] = 4sin 2 (0,) = c 2 UP . (13) 

To obtain the lower bound we shall use the optimizing function u(s) as a a suitable trial 
function in Eq. We obtain 

»l/6 



2 ^ , • 2/ 1 \ 4sinf</»*) cos 3 ( 
c l > Asm 2 (6*) + 



eF n (l) + J^ € F n (u(s))/s 2 ds 



(14) 



e(20* + sin(20,)) 

Since F n is negative, we may combine Eqs. ffTB"]) and f|T4l) and write our main result as 
given in Eq. (J5]). 

As an example consider the reaction term studied by Brunet and Derrida, f(u) = u — u 3 . 
The lower bound can be written explicitly as 

2 • 2/ - n 2(1 - e 2 ) cos 3 (0*) sin^*) 2esin(0 ! ,) f x h 



c >4sin (0*) — — — / cos (f)(s)ds. 

20* + sin (20*) cos(0*)(20* + sin(20*)) h 

The integral has a long analytic expression which we omit here. From the explicit expression 
above it is not difficult to show that the contribution of the two last terms, which arise from 
the nonlinear terms, are of 0( \ lne|~ 3 |. In figures 1 and 2 we show the bounds together with 
the Brunet-Derrida formula as a function of e. The solid lines correspond to the upper and 
lower bounds. The dashed line is the Brunet-Derrida formula. 
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FIG. 1: Speed as a function of the cut-off parameter for the reaction term f(u)=u — u 3 . The 
solid lines correspond to the bounds, the dots to the Brunet-Derrida formula. 
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FIG. 2: As in Fig. 1 for different values of the cut-off parameter. 

As a second example we consider the Fisher reaction term f(u) = u(l — u) with a cut-off. 
The lower bound becomes 

8 '•1ACOS 3 0(S) 



c 2 > 4sm 2 (<f>*) - - 



'1 — e) sin( 



cos" 



'esm 



-d s. 



3 (20, + sin(20*)) 3 (20, + sin(20*)) h 

Again, the integral can be done analytically and we do not show it here. 

In Fig. 3 we show the upper and lower bounds and the Brunet-Derrida formula. In this 
case the Brunet-Derrida formula leaves the allowed band at larger value of e. In general 
for reaction terms f(u) = u — u n , the gap between the upper and lower bounds becomes 
narrower and the Brunet-Derrida formula valid for a smaller range of e. 
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FIG. 3: Speed as a function of the cut-off parameter for the reaction term f(u) = u — u 2 . Lines as 
in Fig. 1 

In summary, we have studied the effect of a cut-off on reaction terms which satisfy the 
KPP condition Eq.(j2J). We have found upper and lower bounds valid for all values of the 
cut-off parameters, which allow to assess the accuracy of the Brunet Derrida formula. If we 
consider only the linear terms, the upper and lower bounds coincide and give the exact linear 
value for the speed, of which the two leading order terms are the Brunet-Derrida formula. 
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